Entanglement properties of optical coherent states under amplitude damping 
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Through concurrence, we characterize the entanglement properties of optical coherent-state qubits 
subject to an amplitude damping channel. We investigate the distillation capabilities of known error 
correcting codes and obtain upper bounds on the entanglement depending on the non-orthogonality 
of the coherent states and the channel damping parameter. This work provides a first, full quan- 
titative analysis of these photon-loss codes which are naturally reminiscent of the standard qubit 
codes against Pauli errors. 

PACS numbers: 42.50.Dv, 03.67.Mn, 03.67.Pp 



I. INTRODUCTION 

Entanglement has been established as a necessary re- 
source for the implementation of many useful quantum 
primitives - teleportation of unknown states key dis- 
tribution ^ and computational speed-up in classically 
exponential- time calculations O H] , to cite but a few ex- 
amples. This valued resource must, however, be pro- 
tected against undesired interactions with the environ- 
ment which lead, ultimately, to the decoherence of the 
quantum system. 

There are two proven ways of safeguarding the fragile 
quantum states from decohering: Quantum Error Cor- 
rection (QEC) codes protect the information through en- 
coding in a larger Hilbert space; Entanglement Purifi- 
cation (EP) protocols aim to distill entanglement from 
a number of identically prepared copies. It is known 
that QEC codes can be recast as EP schemes and vice- 
versa; the connection for discrete-variable (DV) systems 
was established in the seminal paper by Bennett et al 
[5], whereas the corresponding bridge for continuous- 
variables (CV) has been demonstrated in Refs. [BHS]. 

In the optical context, amplitude damping - the ab- 
sorption of transmitted photons - is a dominant source 
of decoherence. It can be appropriately modeled by hav- 
ing the signal interact with a vacuum mode in a beam 
splitter with the appropriately chosen transmissivity pa- 
rameter. Amplitude damping appears thus as a Gaussian 
error, and as such it is known that Gaussian resources - 
the set of operations which can easily be implemented 
through Gaussian ancilla, beam-splitters, phase-shifters 
and homodyne measurements - are of no use in protecting 
the signal state [T TO , and more elaborate (and possibly 
non-deterministic) non-Gaussian operations must be ac- 
counted for. 

In the present paper we review such a QEC scheme as 
proposed by Clancy et a/.^, which protects arbitrary 
coherent-state superpositions (CSS, also known as "cat 
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states" [12 ) through the use of non-Gaussian encoding 
operations. We provide a hitherto inexistent quantita- 
tive analysis of this code's performance, examining its 
entanglement distillation capabilities through Wootters' 
concurrence [T3^. This study is carried both directly and 
through the use of entanglement evolution equations [14 , 
thus demonstrating their prowess over non- Gaussian CV 
carriers - even though the logical information may be 
codified in discrete qubits. 

Entangled coherent states subjected to an amplitude 
damping channel were also investigated by P. Munhoz 
et al. with similar methods [15], however, employing a 
different class of states and aiming towards distinct ap- 
plications from those presented in this contribution [? ]. 

A coherent-state defined as an eigenvalue of the 
annihilation operator {a\a) = a\a))^ can be expressed in 
the Fock (number) basis as 

oo ^ 
n=0 ^ 

A detailed introduction to the properties of coherent 
states can be found for instance in jT6[ [17] . 

Following [TT], we identify the logical qubits as |0)l = 
I — a) and = in the so-called (— ,+) encoding. 
Qubits can equally be defined in the (0, a) encoding as 
|0)l = |0) and = |2a), but it can be shown that 
both are equivalent in their decoherence properties and 
can easily be translated via displacement operations, so 
the first convention will be adopted here. An arbitrary 
qubit superposition is therefore represented as 

\Q^) = ^l=(a\-a)^h\a)) (2) 

where \a\^ + \h\^ = 1 and N{a) is a normalization con- 
stant, N{a) = l + e-2|"l'(a6*+a*6). It is argued that, for 
sufficiently large values of a, | — a) and \a) are approxi- 
mately orthogonal, and N{a) ~ 1; however, present-day 
technologies only achieve limited a sizes ( "Schrodinger 
Kittens" [E]). Therefore, a significant amount of non- 
orthogonality must be considered (fig. fl]). 



2 




FIG. 1. (Color online) Overlap between two coherent states, 
{a\ — a) as a function of the size of the coherent state \a\. 
Inset shows range accessible with current technology. 
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II. AMPLITUDE DAMPING 

Photon loss is considered to be the predominant source 
of errors to affect qubits in the optical context [19 . We 
model such loss by interacting the signal with a vacuum 
mode |0)^ in a beam splitter of transmissivity 77, resulting 
in 



\Qh 



(a| - ay/7])\ - a^J\ - r])i + 



h\oi.Jr])\a^\ - r\)i) 



(3) 



The final state after transmission is obtained by inte- 
grating over the loss mode (denoted here by 

P=lld^P i{m)TT{Q\P)i (4) 

For a single coherent state, this integration is trivial 
and amounts to an amplitude contraction, remaining in 
a pure state. For a superposition, though, the resulting 
state after tracing out the loss mode is now mixed. One 
obtains 

p = (1 - pe)^^={a\ - a^/rj) + b\a^/r])) x H.c. 



1 



= (a| - ay^) - b\a./ri)) x H.c, (5) 



where H.c. is the Hermitian conjugate of the previous 
term and N\a) = 1 - e~^l"l (a6* + a*6). In order to 
simplify the analysis of the decohered state, Eq. ([5| can 
be cast into a more convenient form [11 , namely, 

P = {'^-Pe)\Qa^){Qa^\ ^ PeZ\Q a^) {Q a^\Z (6) 



where Pe = |(1 — e 



-2{l-v)\oc\' 



) is the probability that the 



Pauli Z operator (Z(a|0) l + = <^|0)l — was 

applied (fig. |2|. With this expression, photon loss can 
be seen as having a two-fold effect: first, the amplitude 
of the states is unconditionally reduced from a to ay/fj] 
second, with probability Pe 7 the qubit suffers a phase flip. 



FIG. 2. (Color online) Phase flip probability pe as a function 
of the coherent state size |a|, for channel transmissivities 77 = 
2/3 (blue line, above) and 77 = 0.9 (red line, below). 



[20^ can be used to protect the qubit. Such a code can be 
implemented [11 in the optical setting by sending the in- 
put state through a sequence of three beam-splitters fol- 
lowed by Hadamard gates - a highly non- Gaussian oper- 
ation which implements, up to a normalization constant, 
|0)l ^ |0)l + |1)l, |1)l ^ |0)l - |1)l. The (unnormal- 
ized) encoded state which results is 



-,[\-a) + \a)r^ + h{\-a)-\a)Y 



(7) 



After transmission through the loss channels, another 
Hadamard gate is applied to each of the modes, which 
are then recombined through an inverted sequence of 
beam-splitters. The two ancilla modes are measured 
to provide syndrome information, from which the ap- 
propriate correcting operation can be applied to return 
the signal to its "unflipped" state. Finally, teleport- 
ing the state into an appropriately prepared Bell state 
I — — a) + the amplitude can be restored 

to its original value. 

The three-way redundant encoding achieved by the 
procedure outlined above can correct up to one error; 
therefore, the probability of achieving an error-free trans- 
mission is given by 



P success, 3 = 1 - + 2pg 



(8) 



This can be increased by encoding the input state with a 
higher number of repetitions. We will thus also analyze 
codes with 5, 11 and 51 repetitions. 5-way redundancy 
increases the success probability to 1 — lOpg + 15pg — 6pg; 
a n-repetition achieves 



'-'success^n 



E 



k k 



(9) 



III. ERROR CORRECTION 



IV. ENTANGLEMENT 



Having identified the effect of amplitude damping as 
a phase flip, a traditional 3-mode error-correcting code 



The main focus on the QEC literature cited above lies 
on the achievement of the encodings - for instance, the 
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implementation of the Hadamard gates or the teleporta- 
tion strategy. The scheme's overah performance, though, 
is not quantified except for certain success probabihties 
or the fidehties of the involved operations. However, it 
has been noted that "fidelity is insufficient to quantify 
quantum processes and protocols" [21]. As such, we will 
explore the known fact that QEC codes can be recast 
as EP protocols [5, 22 and, employing the entanglement 
as a figure of merit, provide quantitative benchmarks for 
this codification. 



Direct calculation 



We will consider initial cat states of the form 



IXcKl ,Q!2 ) 
1 



(v^|q^i, 0^2) + e*^Vl - ^1 - Oil, -0^2)) 



with N{ai,a2) = 1 + 2cose^w(l - w)e-'^\'^^\ "^l^^l 
and |ai,a2) a shorthand notation for |ai)|a2)- 

In the case of direct transmission, the first mode is 
kept, while the second is sent through the photon- loss 
channel, resulting in 

Pdirect (1 Pe)\Xa,a^/fj) {Xa,a^/rj\ 

+ PeZ \z , (11) 

with Pg being the phase flip probability for a two-mode 
state, adjusted to preserve the normalization of each com- 
ponent, defined as 



1 _ e^l^r _ g2|a|^(l-r7) ^ g2|a|^(l+r7) 



2(l-e4|«r) 



(12) 



If the sender and receiver make use of the 3-mode rep- 
etition, the entangled pair they will share after encoding, 
transmission and decoding is given by 

Pfinal,3= (1 - 3Pg^ + '2Pe)\Xa,ay^){Xa,ay^\ 

+ (3P,2 - 2P!)Z\Xa,aV^){Xa,aV^\Z, (13) 



with the Z operator taken to act, in both (11) and ilSL 



on the transmitted mode. The use of higher order n- 
repetitions results in 

Pfinal,n— Psuccess,n \Xa,a^/fj) {Xa,ay^\ 

(14) 

~l~(l Psuccess,n) Z\Xa,ay^) {Xa,ay^\Z ' 

with Psuccess,n defined as in ([9|, but depending on Pg 
instead of Pe. 

We adopt Wootters' concurrence to quantify the en- 
tanglement; in the case of a bipartite qubit system, the 
concurrence is given by [13] 



C = max{0, y ^1 — y 



where are the eigenvalues, listed in decreasing order, 
of pp. p is the time-reversed density operator. 



{cTy,! Cr2y,2)P*(^2y,l ^ ^2/, 2) , 



(16) 



where ay^i is the Pauli Y operator in the i-th mode. In 
order to cast the density matrix in the appropriate form, 
we make use of an orthogonal basis {\ua)^\va)} (see for 
instance [23]) such that 



\a) = Pa\Ua) + J^a\Va) 
\- a) = Pa\Ua) - J^a\Va) 



(17) 



with Pa 



1 



-2\a\-' 



and 



-2\a\^ 



(10) This way, the entangled pair can be written, up to a 



normalization factor of 1/ y N{a, a) (omitted below for 
clarity), as 



\Xa,a) = Vui\a)\a) + e^Vl - ^| - a)\ - a) (18) 

= (V^ + e'^Vl - W) {pl\Ua)\Ua) + ^l\Va)\Va)) 
+ {^/w -e'^Vl - w)paJ^a {\Ua)\Va) + \Va)\Ua)) 



and thus 



IXq;,^) (Xq!,q; I — 



(19) 



where 



/ |ap a6* ac* ad* 

6a* |6p 6c* 6d* 

ca* cb* |cp cd* 

Vda* d6* dc* \d\^ 

a {y/w + e'^\/l-w)pl^ 

b = C = {^/w — e*^Vl — w)paJ^a 



After transmission, a similar matrix is also obtained for 
the unflipped state (|Xa,ay^)); the phase-flipped states 
after transmission are described by 



f |ap a6* ac* aJ* ^ 
id* |6P be* id* 



cd* cb* 



cd* 



Z\Xa,ay/r]) {Xa,ay^\Z 

\dd* db* dc* |Jp 
with Z taken to act on mode 2, and 

a (V^ - e'^Vl - w)paPa^ 

i {\/w + e'^Vl - W)paiya^ 

c {^/w + e'^Vl - w)po,^Vcx 

d {^/w - e'^Vl - w)UaJ^a^ . 



(20) 



J 



With (19) and (20), one can construct the matrices 



A4} (15) (11), (13) or (14) and, through the use of (16), calculate 
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the concurrence as defined in ([T5|. Fig. [s] plots tiie en- 
tanglement of the initial state with w = 1/2^ correspond- 
ing to "genuine" Bell states. It is seen that ^ = and 
= 7T (respectively "even" and "odd" due to the parity 
of the number states in the corresponding superposition) 
result in maximum entanglement, the former asymptoti- 
cally for large superpositions, the latter independently of 
the superposition size. 



by an X matrix 




FIG. 3. (Color online) Concurrence for the initial state as 
a function of the size of the coherent state \a\ and the phase 
parameter 0. 

The above method, however, is not convenient for cal- 
culating the entanglement of different levels of encoding, 
as it requires the eigenvalues to be computed for each 
variation, which can be a time-consuming task if done 
analytically. Thus we present, in the following section, 
an alternative way to obtain the state's concurrence. 



B. Entanglement Evolution 

The effect of an arbitrary quantum channel $ to a state 
can also be described in the dual picture [6 , interchang- 
ing the roles of the channel and the initial state. An 
evolution equation is obtained [M], which equates the 
entanglement (and in particular the concurrence) of the 
final state to the product of the concurrence of a maxi- 
mally entangled state subjected to the same channel 
times the concurrence of the initial state |x). The prob- 
lem is thus reduced to the calculation of two, possibly 
simpler, concurrences. 

C [(1 ® $) |X){X|] =C[{1® $) |0+){</.+ |] C [\x)] . (21) 

As in the previous section, the states will be writ- 
ten in the orthogonal basis {\ua)^ \^a)}- One possible 
maximally entangled Bell state in this basis is given as 
\^t,a) = {Wa)\ua) + \va)\va)) / V^, but any statc carry- 
ing exactly one ebit will have its entanglement affected 
in the same way. The resulting state after transmission 
through the amplitude damping channel can be described 



/a /\ 

b z 

z* c 

\f* dj 



(22) 



where 



(1 



.e-2i«r(i-^))^2 



(-l + e-2|«l'(i-^))i/^ 



1.0 

0,5 ConeiuTeiice 
(e1>its) 

0.0 



c = 



_(_l + e-2|a|^(l-.))^2^ 



z = 



Jl. 
(-1 



■e-2|"l^(i-''))^,^^.. 



Vv 



■e-2|"l'(i-'7))/x. 



The outer and inner elements of the (22) represent. 



respectively, "unflipped" and "flipped" Bell states of re- 
duced, ay^, amplitude. Just as in ([5|, this allows one 
to rewrite the resulting state as 

(1 ® $) K J«J = (1 - PeMta^){<l>t^^\ 



Equation \21\ can also be extended to encoded states, 
such as those obtained by the QEC repetition codes de- 
scribed in Sec. III. In this case, the channel $ will in- 
clude not only the lossy transmission channels itself, but 
also the encoding, syndrome measurement, error correc- 
tion and decoding operations. Nevertheless, the resulting 
density matrix for the Bell state is still an X matrix. 



(l®$e„c)KJ«J 



-(1 



success 



^)^<aJZ (?4) 



where P. 



success,n •> 



as in the Sec. Ill, is the probability of 
achieving error-free transmission. 

The concurrence of a state described by an X matrix 
can easily be found by [25J 



(25) 



C {p) = 2 max 0^\z\ — V ad, \ f\ 



The only remaining step is the calculation of the second 
term of the RHS of (21). Through the methods outlined 



in the previous subsection, one finds: 



C 



2(l-e-4|"l') y/w{w - 1) 
1 + 2V'^(^-I)e-4|«I' 



COS( 



(26) 
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This method matches the results found in the previous 
subsection, but ahows for a computationahy significant 
speed-up in calculation times. The concurrences for dif- 
ferent encodings are compared below: 

Concurrence 




(b) r] = 0.9, = 



Concurrence 
1.0 r 




(c) T] = 2/3, 6 = 71 

Concurrence 




(d) r] = 0.9, = 7T 

FIG. 4. (Color online) Concurrence after transmission as a 
function of the coherent state size for different values of 
transmissivity 77 and the phase parameter 0. Direct trans- 
mission (blue, thick line) is compared to encoding with 3 (red 
line), 5 (green, dashed), 11 (black, dotted) and 51 (grey, dash- 
dotted) qubits. 



Two distinct observations can be made from fig. [4] 
First, for the even cat states (a) and (b), the non- 
orthogonality of the basis states prevents the pair from 
achieving high entanglement for low values of a; higher 
encodings are of little advantage in this regime. The 
choice of the phase is of paramount importance when 
one compares to the odd states (c) and (d), where even 
for very small a there is still a significant amount of 
shared entanglement, and encoding in more qubits still 
improves on this amount. Second, for sufficiently large 
sizes of a, the fiip probability approaches 0.5 and thus it 
dephases the qubit entirely, independently of the phase 
or encoding. However, for a certain range, the different 
codifications are seen to help achieve and sustain higher 
entanglement between the shared pairs. 

Finally, in fig. |5]one can observe that, all the way to 
zero transmissivity, the pair has some residual entangle- 
ment. Higher redundancy improves this figure; we also 
note that, in the context of [21 , the encoding does not 
induce entanglement sudden death (ESD) [25] . 




FIG. 5. (Color online) Concurrence after transmission as a 
function of the channel transmissivity 77, shown without en- 
coding (blue, thick line) and encoded with 3 (red line), 5 
(green, dashed), 11 (black, dotted) and 51 (grey, dash-dotted) 
qubits. Here a = 1.3 and ^ = 0. The behaviour for ^ = tt 
follows qualitatively the same form. 



V. DISCUSSION AND CONCLUSIONS 

In this work, we have studied the decoherence pro- 
cess of entangled coherent-state superpositions in the 
amplitude damping channel, exploring QEC codes orig- 
inated from the discrete- variable regime in an optical, 
continuous-variable setting. The quantitative analysis 
was reminiscent of EP protocols, however, this trans- 
lation must be carried out carefully: a true distilla- 
tion would, in addition to the non-Gaussian Hadamard 
gates necessary for the encoding, require teleportations 
to be performed via non-Gaussian two-mode Bell mea- 
surements 13 [24]. In this picture, multiple copies of the 
entangled resource |^+) are shared across the amplitude 
damping channel. The use of the encoding and decod- 
ing circuits of the original error-correcting code then re- 
sults in one pair with higher entanglement (see fig. [6|, 
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which effectively characterizes this scheme as a one-way, 
deterministic entanglement distillation for noisy, non- 
Gaussian CV states. 

Alternatively, in the context of [22^, the translation 
to an entanglement purification protocol could be per- 
formed through an elaborate non-Gaussian multi-mode 
projective measurement in the encoding basis. The non- 
Gaussian elements are, in any case, a required condition, 
enabling the circumvention of known No- Go theorems 
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channels 



decoding 
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(10$)l 
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decoding 



FIG. 6. (Color online) The QEC code as presented in 
(above) and a possible interpretation as an EPF (below). 

The protection granted against phase-fiip (or Z) er- 
rors by different repetition codes was investigated and 
quantified through concurrence; however, possible physi- 
cal limitations in the implementation of the gates neces- 
sary for these encodings have not been been taken into 
account, therefore establishing the above results as an up- 
per bound. In particular, we note that both Hadamard 
and teleportation operations in pT] only approach unit 
fidelity and/or high success probabilities for large super- 
position sizes. However, as we showed here, in this regime 
these codes would cease to work. We conclude that, in 
order to incorporate the actual encoding operations into 
the above protocol, gates that succeed with high fidelity 
on not too-large a cat states will be needed. One possi- 



bility to achieve this could be the scheme of [12], where 
a ~ 1.2 is enough for performing fault-tolerant quantum 
computing. 

We have observed a trade-off between the orthogonal- 
ity of the coherent-state basis and the probability of a 
phase-fiip error. The former imposes a necessity for a 
maximum size of coherent-state superposition, while the 
latter prevents the use of arbitrarily large superpositions, 
thus hinting at an optimal regime depending on channel 
parameters. 

The generation of optical cat states remains an exper- 
imental challenge, though much progress has been re- 
cently achieved. Arbitrarily large squeezed cat states 
can be obtained through photon number states and ho- 
modyne detection [26], however such technique is highly 
probabilistic and presents only moderate output fideli- 
ties. A different approach involves tapping squeezed vac- 
uum in a beam splitter (BS), with one output mode 
directed to a number-resolving photon counter. Condi- 
tional on the number of photons detected, the other mode 
is projected into an odd or even cat state [HI [27]; an n- 
photon subtraction is described as a^6'(eo)|0), where a is 
the annihilation operator, S is the squeezing operator and 
eo amounts to the degree of squeezing applied. Superpo- 
sitions with a ~ 1.2 — 1.3 can be created through this 
technique; it is also worth noting that the more "valu- 
able", 1-ebit odd CSS resource is more readily obtained 
than its less entangled even state counterpart. 

The choice of "odd" cat states (with 6 = n) could be 
seen as a viable alternative to employ small-amplitude su- 
perpositions to distribute entanglement while minimizing 
the amplitude damping effects. Still, a may not be cho- 
sen too small so that the resulting qubits become imprac- 
tical for information transfer or computational purposes. 
Nonetheless, in the context of optical, realistic quantum 
communication, a combination of "odd" CSS resources, 
linear optics and nonlinear measurements for encoding 
may lead to efficient error correction and entanglement 
distillation strategies. 
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